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Abstract. We study the one-dimensional spin- 1/2 antiferromagnetic Heisenberg model exposed to an ex- 
ternal field, which is a superposition of a homogeneous field /13 and a small periodic field of strength hi. 
For the case of a transverse staggered field a gap opens, which scales with hl^ , where ei = ei(/i3) is given 
by the critical exponent rji(M(h-i)) defined through the transverse structure factor of the model at h\ = 0. 
For the case of a longitudinal periodic field with wave vector q = tt /2 and strength hq a plateau is found 
in the magnetization curve at M = 1/4. The difference of the upper- and lower magnetic field scales 
with /i3 — /13 ~ h'^q , where es = e3{hs) is given by the critical exponent ?73(M(/i3)) defined through the 
longitudinal structure factor of the model at hq — 0. 



PACS. 75.10 -b General theory and models of magnetic ordering 



1 Introduction 

The properties of the one-dimensional (ID) spin-1/2 anti- 
ferromagnetic Heisenberg model (AFH) with nearest neigh- 
bour coupling: 



with 



qa{M) = TT 



1 



: a = 1,2 



1 - 2M : a = 3 



(1.7) 



H(/l3) 

Ho 



Ho - 2/1383(0), 

N 

2^ S; • S;+i, 



1=1 
N 



S,(q) = ^e''«5r, a = 1,2, 3, 



(1.1) 
(1.2) 

(1.3) 



Conformal field theory describes the critical behaviour 
at the soft modes In particular the field 

dependence of the //-exponents: 



Va{M) = 



7rv{M) 



(1.8) 



1=1 



in the presence of a homogeneous external field of strength 
hs are well known: 

1. There is no gap. The magnetization curve M = Mlhs) 
is a monotonically increasing convex function 

for /13 > 0; in particular there is no plateau. 

2. In the presence of the field /13 the ground state \ps, S) 
of H(/i3) has total spin S* = 5'j. = NMiJi'i) and mo- 
mentum = 0, TT - depending on S and N . 

3. The low energy excitations which can be reached from 
the ground state \ps, S) by means of the transition op- 
erators S3(g) and S±(q): 

0.3(9, h^) = E{p, +q,S)- E{ps, S), (1.4) 
UJ±{q, hs) = E{ps +q,S±l)- E{ps,S) ± hs, 

(1.5) 



has been computed by means of the Bethe Ansatz so- 
lutions for the energy differences and the spin wave 
velocity |,0 

v{M)^^ lim N[E{ps + 27T/N,S)-E{ps,S)] (1.9) 

4. The 77-exponcnts govern the finite-size behaviour of the 
transition amplitudes: 



{S±l,p, + Tr\S±{7r)\S,Ps) 



N~*oo 



{S,Ps + q3\SM\S,Ps) ^-^iV-3(/.3) (i.u) 



^«:i(/l3) (2,10) 



with 



Kaihs) = 1 - 



(1.12) 



vanish at the soft mode momenta qa = qa{M): 



[2a{M) = lim NuJaiqaiM),h3), 



(1.6) 



and of the static structure factors: 

{S,p,\SMSa{qaV\S,Ps) ""-^ 7V2-'74A^). (1.13) 

At the soft mode momenta qa = qa{M) the dynamical 
structure factors develop infrared singularities of the 
type w-2'^»(''3). 
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First evidence for the existence of low energy modes in 
the excitation spectrum has been found recently in neu- 
tron scattering experiments on copper benzoat, |ll| , |l2|| ex- 
posed to a homogeneous magnetic field /13. An exponential 
fit to the temperature dependence of the specific heat data 
revealed, how e ver, that there is a gap in the energy dif- 
ferences (pr^),(1.5) and (1.7) , which opens with the field 



strength /13 as e = 2/3. This means of course, that 
the compound copper benzoat can not be described by a 
ID Heisenberg antiferromagnet. Oshikawa and Affleck |0| 
argued that the local g-tensor for the Cu ions generates 
an efi^ective staggered field of strength {hi ^ ^,3), perpen- 
dicular to the uniform field /13. Therefore, one is lead to 
investigate the Hamiltonian: 



H(/^3,/»i) =H(/i3) + 2/iiSi(7r) 



(1.14) 



It is the purpose of this paper to study the evolution of 
the gaps 



L0a{qa,h3,hi) CX /ij' 



Xh3 



(1.15) 



by switching on the transverse staggered field. In partic- 
ular we are interested in the /13-dependence of the expo- 
nents Eaihs). 

It has been pointed out by the authors of Ref. |Q 
that a staggered field alone, i.e. = 0, M = 0, generates 
a ground state gap which opens with hi, e = 2/3. In 
a previous paper we have studied the finite-size scaling 
behaviour of the gap and of the staggered magnetization 
in the scaling limit hi ^ 0, N ^ 00 and fixed scaling 
variable x = Nh\ at M = 0. 

The method used in Ref. ^ is based on a closed set of 
differential equations, which describes the /ii-evolution of 
the energy gap uj^jir, 0, hi) [Eq. ( |l.4| )] and of the relevant 
transition amplitude (1.11) f or h^ = 0. It turns out that 
the exponent e(/i3 — 0) in ( [1.15[ ) is fixed by t he fin ite- 
size behaviour of the initial values, i.e. ( |l.4| ) and (1.11) for 
h3^hi= 0: 



e{h3 = 0) 



1 



1 + k(0) 



(1.16) 



to the case 



In this paper, we extend the method of Ref. 

/l3 > 0. 

In section B we discuss the evolution equations for the 
Hamiltonian ([1.141) . The finite-size behavio ur o f the initial 

0) for the gaps ^ and 



> 



conditions {hi — O^h^ ^ ^ ^ 

and fo r the relevant transition matrix elements ( |l.lO[) and 



1.11) is reviewed as well. 



Switching on the transverse staggered field in (1.14) a 
gap opens at the field independent [Eq. (1.5) for q = tt] 
and the field dependent [Eq. ( [T^ for q = 93 (M)] soft 
modes. The finite- size scali ng behaviour of these gaps is 
studied in sections |2.1| and |2.2|, respectively. In section ||, 



we investigate the effect of a longitudinal periodic field on 
the low-energy excitations of the AFH model. From these 
results we infer in section 3.1 the corresponding magneti- 
zation curve. 



2 Evolution equation and initial conditions 

Starti ng fr om the eigenvalue equation of the Hamilto- 
nian (1.14) 



H(/i3, hi)\Wn{h3, hi)) = En{h3, hi)\^n{h3, hi)), (2.1) 

it is straight forward to derive the following set of differ- 
ential equations 



dhl 

dhl 



(2.2) 



dhl ' 
(2.3) 



which describes the evolution of the energy eigenvalues 
En ^ En{h3,hi), energy differences a;„™ = a;„„(/i3, hi) = 
En — Em and transition matrix elements 

Tnm{h3,hi) = {^n{h3,hi)\Si{n)\^.ra{h3,hi)), (2.4) 

of the perturbation operator Si(7r).|^The latter has the fol- 
lowing properties: It changes the momentum by Ap — n 
and the total spin Sj^ by one unit. Therefore, the eigen- 
states |!fVi(ft,3, hi)) are linear combinations 

\^n{h3,hi)) = J2 [an{ST,hi)\p„,S^) 

+bn{Sl.,hi)\pr,+n,S^)], (2.5) 

of eigenstates \pn,Sj^) and \pn + TT, Sj^) to the total spin 
Sj- and the mo ment a p„ , p„ -I- tt. Note, that the evolution 
equations ( ^.2| ) and (2.3) decouple for diff eren t m omenta 
Pn,Pm with \pn — Pm\^ scction 2.1 and 2.2 we will 
study the following cases: 

1. Pn = 0, TT, 

2. p„^q3{M), q3{M)+n. 

For both cases we have the initial conditions: 

(lnm{h3) 



Wnm(<7, h3, hi = 0) 
Tnm{h3, hi = 0) 



N ' 

bnm{h3)N''^'''\ 



(2.6) 
(2.7) 



which are completely fixed by the excitation energies and 
transition amplitudes of the unperturbed problem {hi — 
0) in a uniform field /13. We can now repeat the whole 
line of arguments, we developed for h^ = in Ref. [^. 
The evolution equations (|2.2| ) and (^) possess scaling 
solutions: 



Wnm(g, h3, hi) 



T„m{h3,hi) = Nh'{'''^'^^0nm{x), 



(2.8) 
(2.9) 



^ The A^-dependence of eigenvalues and transition matrix el- 
ements is always understood. 
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in the combined limit 

hi-^0, N oo, 



X EE Nhf'""^ fixed. 



(2.10) 



The exponents £(^.3) and <T{h^) are given by t he finite-size 
behaviour of the initial values (2.6) and (2.7): 



1 



1 + k(/i3)' 



1 - /t(/i3) 
1 + k(/i3)' 



(2.11) 



2.1 The gap at the field independent soft mode g = tt 

As was pointed out in the introduction, the ground state 
\n^{))^ \p,,S) of the ID spin-1/2 AFH model, H(/i3,0), 
in the presence of a uniform field /13 has total spin — 
S = NM{hz) and momentum ps = 0, or ps = tt. The 
first excited state which can be reached with the operator 
Si(7r): 



±l) = \ps + n,S^T^S±l), 



(2.12) 



has a gap of the type (p.6| 



w±io(^, /i3, 0) = E{p, + TT, 5 ± 1) - E{p„ S) T /J3, 

(2.13) 



which vanishes as N 
elements: 



for N 00. The transition matrix 



T±io(/i3,0) = (±1|S±(^)|0) 



N- 



diverge in the limit N — s- 00, where Kii^h^) is obtained 
from the known rji{M) exponent (1.12). Both curves, 771 = 
?7i(M) and M = M{h^) were determined exactly by means 
of Bethe ansatz solutions on large systems § , as well via a 
solution of a system of nonlinear integral equation derived 
from the Bethe Ansatz ||l^. The /i3-dependence is shown 
in Fig. 2.1. It starts at the known value ei{h^ = 0) = 2/3 
and then drops monotonically with h^. At h^{M = 1/4) = 
1.58 ... , the exponent is reduced to 



ei(;i3(l/4)) = 0.5975... 

In order to explore the scaling behaviour ([2 
we have determined numerically the ratios 



(2.15) 
I) of the gap. 



^io(7r, /t3,fel) 
wio(7r, /i3,0) 



= 1 



= 1 



eio{x, hs), 

^io(a;, /13) 
X , 

aw 



with X = Nh'i''^^^'^ and f2iQ as given in Ref. 



(2.16) 



on fi- 



nite systems. The homogeneous field /13 has to be cho- 
sen carefully. According to our premise, the ground state 



|0) = \ps,S) (at hi — 0) has total spin Sj 



S 



NM{h3) and energy E{ps,S) - 2h3S^. The two excited 
states |±1) = \ps+Tr,S^ = S±l) have a gap. Positivity of 
the gaps yields an upper and lower bound of the /13-field 
m >h3> hi): 



2hl 
2hi 



E(ps+Tr,S+l)~E{ps,S), (2.17) 
E{ps,S)-Eips+TT,S-l), (2.18) 



0.047 0.101 0.166 0.254 0.5 




Fig. 2.1. The exact critical exponents ei (solid line) and £3 
(dashed line) versus /13 and M, determined from a Bethe ansatz 
solution of finite system size A'' = 4096. 



which leads to the well known steps in the magnetization 
curve on finite systems 0. Note, that at the edges /if and 



/ig the excitations energies: 



LJ+l0i7T,h^,hl^0) = 0, 

w_io(7r, = 0) = 0, 



(2.19) 
(2.20) 



vanish identically. Therefore, ratios of the gap (2.16) do 
not make sense in these cases. At the midpoint field /13, 
however: 

2h3 = {hi + h\)/2 

= [E{p,+^,S+l)-E{ps + T:,S-l)]/2, (2.21) 

the two excited states have the same gap: 



w±io(7r, /i3,0) = 

E{p,+^,S+l) 



E{p,+Tr,S-l)-2E{ps,S) 



(2.22) 

The degeneracy of these two excited states is not lifted 
in the first oder perturbation theory in hi, since all the 
relevant matrix elements 



(n|Si(7r)|TO) = 0, n,m = ±l 



(2.23) 



vanish. The ratio ( |2.16| ) is shown in Fig. |]^(a), for the 
midpoint field /13 = /i3(iV) « 1.58, corresponding to a 
magnetization M = 1/4 on system sizes N — 8, 12, 16, 20. 
Optimal scaling is achieved here, with the exponent ei = 
0.595( 5), w hich is in excellent agreement with the exact 
value ( 2.15 ). According to Ref. |9), the low x-behaviour of 
the scaling functi on e in(x, ^3) i s also predicted by th e evo - 
lution equations (2^) and (2_^) in the scaling limit ( ^.10 ): 



eio(a:,/i3) = eio(/i3)^"^^^''^\ 



(2.24) 



with ipiih^) — 2/ei(/i3). The linear behaviour in the vari- 
able for small x- values is clearly seen in Fig. 2.2(a). 
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1.3 



1.2 



1.1 



(a) E|=0.595 



1.4 



1.3 



1.2 



1.1 



1 * 



(b) E,=2/3 



1 2 

2/e, 



Fig. 2.2. A comparison of the ratio ( p.lq ) for two different 
values of ei and the midpoint field fts [Eq. (2.21)] for system 
sizes iV = 8, 12, 16, 20 



The effect of the homogeneous h^-Held on the exponent 
ei is demonstrated in Fig. |2.2| (a). An exponent ei{h^) — 
ei(/i3 = 0) = 2/3 independent of hs w ould lead to consid- 
erable scaling viola tions of the ratios (2.16), as is demon- 
strated in Fig. |^(b). 



2.2 The gap at the field dependent soft mode q = qs 



Let us now turn to the field dependent soft mode [Eq. (1.4) 
for q = (A/)]. Switching on the perturbation operator 
hiSiin) the ground state energy E{ps, S) and the energy 
E{ps+q3{M), S) of the excited state evolve independently, 
since their momentum difference (?3(M) does not fit to the 
momentum transfer tt mediated by the operator Si(7r). In 
other words, we have to study the ground state energy 
£'0(^3? ^1) in the sectors with momentum p = 0, tt and p = 
q3{M), q3(M)+ Tr, s eparately. In both cases insertion of the 
scaling ansatz (|2.8|) and {2S) for the e xcita tion energies 
LOno and transition amplitudes r„o into ( |2.2| ) yields: 



^^^0 ^ _^1+2ki(/i3)2.1-2ki(/!.3) 

dh'i ^ 

1 



\em{x)? 
f2iQ{x) 



(2.25) 



where x = Nh\^''^''\ To integrate ( ^!25| ) we introduce 

y ^ ^iAi(h3) ^ ;j^^i+«i(/i3)^ (2.26) 



and 



(2.27) 



from which follows: 



Eoihs, hi) 



Eo{h3,0) 

ei(/i3) rV 



dy' / dy"!{v"). (2.28) 



Here, we have used the fact, that 



dEo{h3,hi) 



dhi 



= (0|Si(^)|0)|,.,=o = 0. 



(2.29) 



hi=Q 



Equation ( 2.28 ) describes the lowering of the ground state 
energy, if we switch on the staggered field of strength 
hi. We observe the same scaling behaviour with h^^ , we 
found for the excitation energies w„o(7'') hs, hi). In Fig. 2.3 



we have plotted [i?o(^3, ^1) - £^0(^3, 0)]/hl^^'^^' versus the 
scaling variable x-^^'^^~^ for the case p = 0, tt. We observe a 




Fig. 2.3. The scaling of the g round state ene rgy ( | 2.29| ) for the 
midpoint field hi [Eq. (|^)] in the limit ( |2.10 ). Numerical 
data were obtained on system sizes N — 8, 12, 16, 20. 



linear behaviour in this variable, which is a consequence of 
the small x-behaviour of the energy d iffere nces SIiq (x) and 
transition amplitudes 6';o(a;) in Eq. (2.27) [See Ref. M]: 



(2.30) 



Therefore, the integrand on the right-ha nd sid e of ( 2.28| ) is 
constant and the small x-behaviour of ( 2.28 ) is governed 
by y2-ei = x^/'^i-i. 

Let us next turn to the lowering of the ground state 
energy in the sector withp = q3{M), q3{M)+TT. The expo- 
nents K±(/i3) are defined by the initial conditions (hi — 0) 
for the transition matrix elements: 

(±llSi(7r)|0) = {p,±i+q3{M),S±l\Si{7r)\ps+q3{M),S) 
= b±io{h3)N^^<^'''\ (2.31) 

Confornial field theory relates the corresponding yy-expo- 
nents {k± 1 — ri±/2) to the scaled energy differences: 



V±{M) 



n±{M) 

nv{M) ' 



(2.32) 



with 



n±{M)^ lim N[E{ps±i+q3{M),S±l) 

N^OQ 

-E{ps+q3{M),S)]. 



(2.33) 
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Here v{M) is the sp in wa ve vel ocity (1.9) at the soft mode 
q = 0. Evaluating ( 2.33 ) and (L9) l eads to the following 



representation of the ?7± -exponents ( |2.32| ) 

v±iM) 



77±(M)=77i(Af) 



v{M) 



(2.34) 



turn out to be finite, except for the case, where we meet 
a soft mode: 

UJ3{q, /l3, ^ = 0) = E{ps+q, S, hq ^ 0) ^ E{p,,S, hq = 0) 
JV— oo a3(/l3) 



N 



where 



This happens if: 



v±{M) = — lim N[E{ps±i+q:i{M±l/N)±2TT/N,S±l) 
-E{ps±i + q3{M±l/N), S ± 1)], (2.35) 

are the right-hand- (+) and left-hand (-) spin wave ve- 
locities obtained from the slopes of the dispersion curve 
approaching the soft mode momentum from the right- and 
from the left-hand side, respectively: 



g = g3(M) = 7r(l-2M), 



(3.4) 



(3.5) 



e.g. a soft mod e ap pears at g = ±7r/2 if M = 1/4. At 
the soft mode (3.5) the transition matrix elements ( p^ ) 
diverge: 



T3(/l3,0)^^63(/l3)^"^'"^' 



(3.6) 



p^Ps+q3{M)±27:/N. 



(2.36) 



From confo rmal invariance arguments for the energy dif- 
ferences in (2.35) we get 



with an exponent K3(/i3) = 1 — ry3(A/(/i3))/2, given by the 
773(M)-exponent, given in the introduction. Fro m the evo- 
lution equations with the initial conditions (3^) and ( |3.6|) , 
we get in this case a finite-size scaling behaviour of the gap 
ratio: 



r)+{M) ^r]^{M) = l + r/i(M). 



(2.37) 



In summary, we conclude that the gap of the field depen- 
dent soft mode q3{M): 

E{ps+q3{M), S) - E{ps, S) ^ hl'^'''\ (2.38) 

is dominated by the lowering of ground state energy E{ps, S) 
and therefore scales with the same exponent eiihs) as the 
field independent one. 



3 Opening of a gap in a longitudinal periodic 
field 



^3(93, ^3, hq) 
^^3(93, ^3,0) 



1 + 63(2;, ^,3), 



(3.7) 



with a scaling variable x = Nhq^^'^^\ where e^i^a) — 
l/[l-t-K3(/i3)]. The curve e3(/i3) is shown in Fig. 2A. Note, 
that ri3{M) = l/rji^M), which means £3(0) — for 
M = 1/4, we have 



e^ihsiM = 1/4)) = 0.81011 



(3.8) 



A test of the finite-size scaling behaviour (3.7) for q = 
7r/2 and M = 1/4 with the exponent (3. J) is shown in 
Fig. 3.1. The small x-behaviour of the gap ratio is properly 
reproduced with a;^/*^^ and compared with th e pr ediction 
h^q^, where £3 = esihsiM = 1/4)) is given by (|]|). 



So far we have only considered the Hamiltonian (1.14) 
with an inhomogeneous field ft,iSi(7r) transverse to the ho- 
mogeneous field ^.38 3(0). By means of the evolution equa- 
tions (|]^) and ( |2.3| ) we can also study the influence of a 
longitudinal periodic field 



U{h3, hq) = Ho - 2/1383(0) + 2hqS3{q). 



(3.1) 



The perturbation operator 83(9) = [83(9) -I- 83(— g)]/2 
commutes with the total spin operator 5^ and changes 
the ground state momentum ps by ±g. For this reason, all 
momentum states with 



1.04 



1.03 



1.02 



1.01 



Pk^Ps±kq, fc = 0, ±1,±2,. 



(3.2) 











□ 










■ 















■ 



□ 





























are coupled via the evolution equation. For example for 
q — 7t/2 there are 4 different momentum states withp^/Tr = 
±1/2, 0, 1, which have to be taken into account. In general. 



2/e, 



the transition matrix elements at /i„ 



0: 



Fig. 

8,12, 



3.1. Finite-size scaling of the gap ratio (3.7), for A'^ 
16, 20 and qsiM = 1/4) = tt/2 with £3 = 0.81 . . . 



n{h3,hq^0) = {ps±q,S\S3i±q)\Ps,S) 



(3.3) 
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3.1 The magnetization curve in a periodic field 

Let us fina lly discuss the influence of the periodic pertur- 
bation in on the magnetization curve M = M{h^). 
First of all one should notice, that the opening of a gap 
for /iq > in the energy differences (3^) does not imply a 
priori a plateau in the magnetization curve. The criterium 
of a plateau with an up per a nd lo wer critical field h 
can be read from ( |2.17| ) and ( ^.1^ ): 



3 J "3 



2{hl-h\)^ lim [E{p,+'K,S+l,hq)-2E{p,,S,hg) 

N~>-oo 



E{p,+Ti,S-l,hq)] 



(3.9) 



The emergence of th e pl ateaus in the magnetization curve 
can be seen in Fig. |3.2| . A fi nite -size analysis shows that 
a non vanishing difference ( |3.9D remains in the thermo- 
dynamical limit. For this analysis we have used the BST- 
Algorithm The /ig-dependence of the plateau width 

is plotted in Fig. 3.3, together with the predicted scaling 
behaviour h^^ for q ~ 7r/2. 




Fig. 3.2. The evolution of a plateau in the magnetization curve 
at M = 1/4, induced by an external field ( [3.l| ) with period q = 
7r/2. The magnetization curve is calculated from finite system 
(A'^ = 20) via midpoint magnetization in conjunction with a 
finite-size extrapolation of the plateau width from system sizes 
of iV = 8, 12, 16, 20. 




Fig. 3.3. The evolution of the difference (3.9) between the 
upper and lower critical field at the plateau M — 1/4. The 
solid line shows a fit to the data for small values of the external 
periodic field hq. Th e ex pected behaviour is oc hq^, with £3 = 
0.8101 given by Eq. (3.8). The dashed line represents a linear 
fit for larger values of hq. 

4 Conclusions 



This paper is aimed to study the effect of a small periodic 
field on the eigenvalue spectrum of the ID spin-1/2 AFH 
model. We are interested in particular in the opening of a 
gap in those situations, where the unperturbed model is 
known to be critical. The critical exponents 771 (M), 773 (M), 
which govern the di verge nce in the transition matrix el- 
ements ( 1.10 ) and ( 1.11 ) of the unperturbed model are 
known. Following conformal field theory, they are related 
to the finite-size behaviour ( |1.6| ) of certain energy differ- 
ences (1.4) and (1.5), which can be computed on very large 
systems by means of Bethe ansatz. 

The evolution of the eigenvalue spectrum under the 
infiuence of perturbation of stre ngth hq is des cribed by a 
system of differential equations (|2.2|) and ( |2.3| ) , which has 



been shown t o hav e scaling solutions (2.8) and (2.9) in the 



scaling limit ( 2.10| ). The exponents e and a in the scaling 



solutions are uniquely determined by the corresponding 
77-exponents in the unperturbed model. We have studied 
in detail the following types of perturbations. 

1. A transverse staggered field together with a homoge- 
neous longitudin al fi eld /iiSjJtt) + /i3S3(0). Both en- 
ergy differences ( |1.4D and (1.5) at the soft mode mo- 
menta (^]^) were shown to evolve a gap with an expo- 
nent 



4-?7a(M(/l3)) 



(4.1) 



with a — 1 depending on the external homogeneous 
field /13 with magnetization M(/i3). 
A longitudinal homogeneous and periodic field 2/1383(0) 
-|- 2hqS3{q). Such a perturbation creates a plateau in 
the magnetization curve M — M{h3) at 



M=ifl-i 

2 V TT 



(4.2) 
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In other words q has to meet the soft mode momen- 
tum q = qsiM) = 7r(l - 2M). The difference of the 
upper and lower critical field, which defines the width 
of the plateau, evolves with an exponent cs^h^), which 
is related to the corresponding rya-exponent via (4.1) 
for a = 3. 
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